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ANALYTICAL STUDY OF PRESSURE BALANCING IN GAS FILM SEALS 

by John Zuk 
Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio 

ABSTRACT 

Proper pressure balancing of gas film seals requires knowledge of the 
pressure profile load factor (load factor) values for a given set of design 
conditions. In this study, the load factor is investigated for subsonic and 
choked flow conditions, laminar and turbulent flows, and various seal en- 
trance conditions. Both parallel sealing surfaces and surfaces with small 
linear deformation were investigated. The load factor for subsonic flow 
depends strongly on pressure ratio; under choked flow conditions, how- 
ever, the load factor is found to depend more strongly on film thickness 
and flow entrance conditions rather than pressure ratio. The importance 
of generating hydrodynamic forces to keep the seal balanced under severe 
and multipoint operation is also discussed. 

NOMENCLATURE 

A area, in.^ 

C constant in friction factor - Reynolds number relation 
C L velocity entrance loss coefficient 

D hydraulic diameter, 2h, in. 

F force, lbf 

F pressure profile load factor (load factor) 

f mean Fanning friction factor 
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h film thickness (gap), in. 

L sealing face (dam) flow length, in. 

M Mach number 

n exponent in friction factor - Reynolds number relation 

P pressure, psia 

AP sealed pressure differential, psi 

Re leakage flow Reynolds number 

T(x) linear tilt factor, h^h^ + ox)/2h m (hj + qx)^ 

u mean velocity in x-direction, ft/sec 

W flow width, in. 

x coordinate in pressure gradient direction 

a relative inclination angle of surfaces, m rad 

/3 film thickness ratio, hj/hg 

y specific-heat ratio 

rj geometric balance ratio or modulus 
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ju absolute or dynamic viscosity, Ibf-sec/in. 

o , 4 

p density, lbf-sec /in. 

Subscripts 

0 sealed (reservoir) conditions 

1 seal entrance 

2 seal exit 

3 ambient sump conditions 

f friction 

i inertia 
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hs hydrostatic 

sd sealing face (dam) 

n net 

m mean 

x flow length location 

INTRODUCTION 

One of the prime objectives in gas film face seal design is to insure 
that the face loading is sufficiently low so high heat generation and high 
wear are prevented; however, contact or close clearance operation must 
be maintained at all operating conditions. Seal balance can be achieved, 
at least theoretically, by properly adjusting the secondary seal diameter 
(see Fig. 1). A common term used by seal designers is the geometric 
balance ratio or modulus. This modulus is defined as the ratio of the hy- 
drostatic closing area to primary sealing face (dam) area and is used to 
determine the location of the secondary seal diameter. It is desirable to 
predict this location analytically. 

Unfortunately, a gas film seal will usually only be ’ 'balanced” at one 
combination of operating conditions. This is the case because the pressure 
profile factor varies with the sealed gas pressure differential. The pres- 
sure profile load factor is defined as the ratio of sealed hydrostatic pres- 
sure closing force to sealing face pressure opening or separating force. 

The pressure profile factor is defined as the ratio of the net or average 
sealing (dam) face pressure to sealed pressure differential. From hereon, 
the pressure profile load factor will be referred to as the load factor. 

Both the load factor and geometric balance ratio have other names in the 
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literature (2) through (7) and sometimes defined in slightly different ways. 
Some of these names are listed in Table 1. Since the load factor can equal 
the geometric balance modulus at only one set of operating conditions, it is 
therefore impossible to completely balance an ordinary face seal for all 
situations. Engineering judgment must be employed to select the proper 
design. Furthermore, since gases are poor boundary lubricants, contact- 
ing gas film seals can generally only carry 1/4-lb to 3/4-lb per inch of 
circumference (1). Since the load factor strongly depends on the sealed 
pressure ratio for subsonic flow, a slight unbalance at large pressure dif- 
ferentials will greatly overload the seal. 

Since advanced rotating machinery may require large pressure dif- 
ferential seals, the seal flow may be choked. Hence the emphasis of this 
study will be on how the load factor varies under choked flow conditions. 
Both parallel and linearly tilted^surfaces and entrance-effects will be con- 
sidered. 

ANALYSIS 

Seal Force Balance Fundamentals 

The importance of the load factor and the geometric balance ratio can 
be illustrated by considering a face seal force balance. The basic equation 
defining seal closing force is (see Fig. 2) 


NET CLOSING FORCE = F g + (F f + Fj) + A RS 


AP - W 



P dx 


[ 1 ] 


Design philosophies differ; however, a common pressure balancing practice 
for gas film seals is to select the spring force F 0 , to overcome only the 

S 7 
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frictional forces F^ and the inertial forces Fj. (The frictional forces are 
due to the secondary seals (e.g. , O- rings, piston rings) and the antirota- 
tion lugs (e.g. , torque pins) rubbing on the housing. ) 

A fundamental consideration in designing pressure balanced seals is 
the selection of the secondary seal diameter. This diameter determines 
the hydrostatic (pneumatic) closing force as illustrated in Fig. 1 . By 
proper positioning of the secondary seal diameter, this closing force can 
be equal to the sealing dam pressure opening force or at least theoreti- 
cally, any degree of seal face loading. The secondary seal diameter can 
be found from the geometric balance ratio, 77 where 


geometric balance ratio, 77 = 


A HS 

a sd 


[ 2 ] 


Another important parameter is the pressure profile load factor, F 
which is defined as the pressure (pneumatic) opening force normalized to 
the sealed pressure differential force acting over the entire sealing dam 
(seal face) area or 


p _ Pressure opening force 
AP A SD 


[ 3 ] 


Across a sealing dam 


Opening force = W / 

JO 


P dx 
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If the seal opening force is equated to the hydrostatic closing force 


W 



P dx = AP A rs 


and substituting this condition into Eq. [3], the result is 


F = 


A HS 

a sd 


= V 


[ 4 ] 


[5] 


When this situation exists, that is the load factor is equal to the geometric 
balance ratio, the seal is said to be perfectly balanced . 

Once the load factor is known, the seal balance diameter can be simply 
calculated. 

Seal Balance Diameter = F(2Rg - 2R^) + 2R^ 


(For a perfectly balanced seal, the seal balance diameter equals the sec- 
ondary seal diameter. ) For some cases the sealing dam opening force can 
be evaluated analytically and hence the load factor can be predicted analyti- 
cally. Examples will now follow. 

Incompressible Fluids 

The simplest case that the load factor can be analytically predicted is 
for laminar, incompressible flow. For purposes of illustration, since 
many seal users are sealing liquids, this case will be developed. The 
pressure profile can be found from (7) 


P = 


p i- 


A Px 


T(x) 


L 


[ 6 ] 



h«(2h 1 + ax) 

Where T(x) = — - - r — , the tilt factor 

2h m( h l + 0!x ) 2 •' 

The net seal opening (separating) force is found by integrating the 
pressure profile across the sealing face. 


where 


F n = WL AP 


(A) 



(See Fig. 3 for an illustration of the symbol usuage. ) 
Thus, the load factor is 



A SD A p 3 + 1 


Where 

Note for parallel surfaces, 0 = 1 

F = 0.50 

(b) Converging faces, 0 = 2 

F = 0.67 


(c) Diverging faces, 0 = 1/2 
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The practical range of load factors appears to vary from 0.6 to 0.9 
for liquid face seals currently used. 

Gas Film Seals 

A. Parallel Sealing Surfaces 

For gases under conditions of compressible viscous subsonic flow, 
the pressure profile for parallel films is expressed as (8). 



By algebraic manipulation, it can be shown that the load factor, 


F = 


n 


A SD AP 3 


1 + 


1 + 


[ 12 ] 


(This formula is identical to the balance modulus presented in (9) and 
the flow coefficient derived for shaft riding seals found in (5). ) 
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Note that for pressure ratios close to one the gas should behave similar 
to an incompressible fluid. This is the case since the load factor is 0. 5. 

At very large pressure ratios 


F = 


lim 

V P 1 - 0 





2 

3 


[13] 


This high pressure limiting situation is not physically correct, how- 
ever. In reality, the governing equations break down in that fluid inertia 
forces become important and the phenomenon of fluid choking can occur. 
This will be discussed in a section to follow. 

Equation [12] is shown plotted in Fig. 4 and illustrates the strong 
variation of the load factor with pressure ratio. Since F = F^g/P^), F 
can equal the geometric balance ratio. only at one operating condition. This 
illustrates why it is impossible to completely balance a face seal for all 
sealed pressure situations and that a slight unbalance at large pressure 
differentials will greatly overload the seal. 

Since Eq. [lOjalso describes the pressure profile for laminar and fully 
developed turbulent flows, the load factor predicted by Eq. [13] applies to 
turbulent flow as well. Note Eq. [13] is independent of film thickness and 
fluid properties. 

B. Sealing Surfaces Represented by Small Linear Tilts 

The pressure profile for sealing surfaces undergoing small linear 
tilts and operating in the compressible viscous flow regime is (8) 
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p = 


1 + 



- 1 


h 2 x(2h 1 


+ ox) 


2Lh m( h l 


+ OX 


[14] 


The load factor can be found by integrating Eq. [14]. The result is 
after letting 0 = h^/hg, and X = Pg/P^. 

For (3 > 1 


F = 


0 


\fl 


0 - 1 


(l - x)(0 - i)^fe 


cos 


-l( \ + 0A 
0 + A 


[15] 


For 0 < 1, two cases are considered since this change in the function 
avoids imaginary terms. 

r 

Case I 


0 2 - X 2 > 0: 



[ 17 ] 
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Note when 0 = 1, the result is the parallel film case which is described 
by Eq. [13]. Eqs. [15], [16] and [17] are presented in (9). (However, in 
(9) the argument of the arccosine in Eq. [15] has a pressure ratio omitted 
in the numerator. The plotted results are correct, however, in (9). ) 

Choked Flow Conditions 

For subsonic viscous flows the analyses yield the analytical solutions 
presented in the previous section. For large pressure ratios and/or 
relatively small flow length-to-gap ratios, the inertia forces must be 
accounted for and choked flow will result. Choked flow in seals generally 
occurs when the physically limiting condition of sonic velocity is reached 
at the seal exit. 

Since inclusion of fluid inertia makes the flow equations nonlinear, the 
flow is solved using an approximate integrated average model. Also, under 
choked flow conditions, sufficiently high velocities occur "at the entrance 
that entrance pressure losses are significant. 

A. Parallel Sealing Surfaces 

The analysis is described in detail in (7). It is separated into two 
parts, which are considered separately and then matched. One part is an 
analysis of the entrance flow, while the other part is an analysis of the 
seal leakage path itself. The governing conservation equations are com- 
bined through relations to yield a Mach number - mean friction factor 
relation. The mean friction factor is related to the Reynolds number by 
the relation 


f = CR‘ n 

C 


[ 18 ] 
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Where the values for the coefficient C and exponent n depend on the flow 
regime as discussed in (8). Once the Mach number at any flow length dis- 
tance x;, is found from the relation 


4fx _ 

1 

1 - “? 1 - Mfl 

> + y + 1 . < 

r y +x m\ y + 1 m^ "i 

In 2 - - in y V 

D 

r 

L M i >4 J 

2y 

L 2 1 2 X JJ 


[19] 


The pressure distribution can be found from 


P x = 



Then the load factor is found from 


/’ 


P x dx 


a sd ap 


[ 20 ] 


[ 21 ] 


The integral in Eq. [21] is evaluated numerically in (7) by using Simpson’s 
rule. 

B. Sealing Surfaces Represented by Small Linear Tilts 

When the sealing surfaces are deformed, the area change complicates 
the equation relating the mean friction factor and Mach number squared. 
The equation that must be solved is 
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dM 2 

dx 


2M11 + 


Y- 1 
2 



h(l - M 2 ) 


fyM 2 f (X) - C h + xsinc, r 


[ 22 ] 


This equation must be solved numerically. The Runge-Kutta solution tech- 
nique used is described in (10). Again, once the Mach number distribution 
is known, the pressure distribution can be found from 



The load factor is again evaluated in (10) by evaluating Eq. [21] numeri- 
cally using Simpson’s rule. 

RESULTS AND DISCUSSION 

Results will be obtained utilizing the various solutions representative 
of both subsonic and choked flow conditions. The first case has already 
been considered (i.e. , subsonic flow between parallel sealing surfaces). 

As previously described, Fig. 4 shows the load factor variation with pres- 
sure ratio varying from one through 10. Figure 5 extends the load factor 
variation to pressure ratios of 1000 for parallel sealing surfaces. Both 
isentropic entrance conditions and conditions representative of an entrance 
loss coefficient of 0.6 are considered. Note for the isentropic entrance 
condition case, the flow is choked for a pressure ratio of ten and the flow 
is probably turbulent when the pressure ratio has exceeded 50. This case 
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illustrates that a maximum limiting value of about 70. 5 is reached for 
ultra high pressure ratios. The results in Fig. 5 were obtained with the 
film thickness fixed at 0. 0002 in. This film thickness is found to be rep- 
resentative of gas film seal operation (4). Also shown plotted in Fig. 5 is 
Eq. [12], which is the classical subsonic flow equation which, although not 
applicable in a strict sense, shows that a limiting value of 0.67 is pre- 
dicted. However, if a 0.6 entrance loss is accounted for (which may be 
the case in real seal applications (8)), Fig. 5 shows that a maximum value 
(0.645) is achieved for choked flow, then the load factor actually decreases 
slightly with pressure over the entire range studied. Thus, as anticipated, 
since the fluid physics is different for choked flow, the behavior of the load 
factor differs greatly than for subsonic flow conditions. 

Figure 6 presents the load factor variation with pressure ratio for 
both linearly diverging and converging sealing faces. Beta ^hj/hg) is the 
parameter that is varied in this study as the pressure ratio is varied from 
one-tenth to ten. Beta equal to one represents the parallel sealing surface 
reference case. This figure also shows that the numerical solution (10) 
agrees with the analytical solutions Eqs. [15], [16], and [17], This figure 
illustrates the deviation of the load factor with pressure ratio which can be 
applied to evaluate the magnitude of the induced lifting force or contacting 
forces due to coning of the seal face. As expected, the larger the defor- 
mation, the more the deviation from the parallel sealing surface case. If 
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thermal coning should be present in a gas film seal, Fig. 6 can be used to 
predict the resulting increase in seal face loading which would be detri- 
mental to the seal operation. On the contrary, if the application permits 
using a converging tapered face, this figure can be used to estimate the 
degree of added or increased separating force that may be generated. 

For subsonic flow conditions, it was seen that the load factor is in- 
dependent of mean film thickness. Under choked flow conditions, the load 
factor does depend on film thickness but not strongly over a practical 
range of film thicknesses. This is illustrated in Fig. 7. The film thick- 
ness is varied from 50 to 500 microinches. A family of curves is shown 
for sealed pressure ratios of 10, 50, 100, 1000, and 10000. The entrance 
flow conditions were assumed to be isentropic. 

Figure 8 essentially repeats Fig. 7 except the entrance flow was as- 
sumed to behave as a flow with a 0.6 entrance loss coefficient. By com- 
paring Fig. 8 with Fig. 7, it is seen that this case more strongly depends 
on the mean film thickness. This dependence, of course, is especially im- 
portant for noncontacting operation which will usually be the mode of op- 
eration for high pressure differential sealing. 

Figure 9 shows the load factor-film thickness variation for choked 
flow and a one milliradian tilted seal face (divergent). For high pressure 
operation the tilt angle is a more meaningful parameter rather than the 
inlet-to-exit film thickness ratio, j3. It can be shown that 


a = 


^m* 1 - f> 

AR(1 + 0) 
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Alpha values of one and two milliradians are representative of face defor- 
mations that can occur in gas film seals. Generally alpha values are so 
small that they are virtually impossible to measure but can be found from 
stress analyses of the seal ring. Figure 9 also illustrates the various 
flow regimes that the seal can encounter. 

Figure 10 shows the corresponding negative one-milliradian face de- 
formation case (converging face). Note the lesser sensitivity to film 
thickness than the divergent face case. 

The load factor variation for relatively large tilted surfaces of posi- 
tive and negative two milliradian and a pressure ratio of ten is shown in 
Fig. 11. Note the extreme sensitivity of load factor on the mean film 
thickness, indicating an extremely high face loading if a means of achieving 
positive surface separation is not provided. 

Pressure Balancing With Hydrodynamic Forces Present 

For gas film seals operating under severe conditions, e.g. , high 
pressure ratios, operation under conditions with rubbing contact may be 
virtually intolerable due to high face loading. Thus, hydrodynamic lift 
devices such as self-acting lift pads (11) must be incorporated into the pri- 
mary seal face design. The seal force balance should now take into account 
this additional force. 

In some applications such as aircraft gas turbine seals, the sealed 
pressure ratio depends on such factors as the mode of operation, e.g. , 
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takeoff, climb, cruise and ambient pressure which varies with altitude. 
Unfortunately, the geometric balance ratio, 77 , is of necessity fixed for a 
particular design. Thus, a suitable secondary seal diameter must be 
chosen which will enable satisfactory operation over a wide range of con- 
ditions. 

Table 2 shows a comparison between the geometric balance ratio and 
the load factor for the NASA self-acting lift pad seal. (The operating film 
thickness effect on the load factor is accounted for. ) The geometric 
balance ratio is, of course, constant for all design points. However, the 
load factor varies (see table 2). And this is because of the different pres- 
sure profile shapes for each design point. Note that the net closing force 
due to the hydrostatic pressure is small (4. 7 to 2.0 lbf) at the four design 
points. Thus, the self-acting pads act, principally, against the spring 
force in this particular design. This is the design philosophy used for the 
NASA self-acting lift pad design. The philosophy could vary with applica- 
tion. 

As previously mentioned in a conventional seal, the net closing force 
is resisted by solid-surface rubbing contact; thus, a total force balance is 
achieved. But in self-acting seals the force balance is achieved without 
rubbing contact. Therefore, for a given design point, the seal will operate 
at a film thickness such that the total opening force exactly balances the 
total closing force. This operating film thickness is obtained by plotting 
total opening forces and total closing forces as a function of film thickness. 
The intersection (see Fig. 12) of these curves is the steady state equilib- 
rium (operating) film thickness. 
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The incorporation of self-acting lift pads can be viewed as a means of 
noncontact operation of a gas film seal under a wide range of conditions. 
These conditions would yield face loadings that are too high for ordinary 
gas film seals. The conditions could be due to pressure, face distortion 
or entrance conditions. Figure 5 shows that entrance conditions can have 
a marked effect on the load factor. Since entrance conditions can vary for 
given designs, the lift pad forces can be designed to accommodate for 
these variations. Thus for high pressure gas sealing applications, incor- 
poration of lift pads are not only important for high axial film stiffness for 
dynamic tracking, but also for maintaining positive film separation. 

CONCLUDING REMARKS 

An analytical study of pressure balancing in gas film seals was con- 
ducted. Both subsonic and choked flow conditions, parallel and deformed 
seal surfaces represented by small linear tilts were investigated. Ana- 
lytical predictors for the load factor in the compressible viscous, subsonic 
flow regime are presented. Conditions of seal operation under choked flow 
were investigated by utilizing numerical solutions. The following pertinent 
results were found. 

1. For classical viscous subsonic flow, the load factor is a strong 
function of sealed pressure ratio and increases with pressure ratio but is in- 
dependent of film thickness. This is true for both laminar and fully de- 
veloped turbulent flow. 

2. For choked flow the load factor depends not only on pressure ratio, 
but also on film thickness, entrance conditions, and flow regime. 
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a. For isentropic entrance conditions, the load factor reaches a 
limiting value of approximately 0. 70 at pressure ratios greater than 50. 

The flow is turbulent for these conditions. 

b. For seal flows with a 0.6 entrance loss coefficient, a maximum 
of approximately 0.64 was reached; then the value began to drop fairly 
rapidly until the pressure ratio exceeded 40 where turbulence was set to 
occur. The load factor then decreased relatively slowly. 

3. The film thickness was varied from 0. 00005 in. to 0. 0005 in. for 
conditions representative of gas film seal operation. Results indicate that 
for pressure ratios greater than 5, the load factor was sensitive to film 
thickness for a fixed pressure ratio. Also the same study with 0. 6 en- 
trance loss conditions showed film thickness sensitivity but differed greatly 
from the isentropic entrance condition case. 

4. The effects of linear tilts representing convergent and divergently 
deformed faces indicated the importance of controlling and accounting for 
seal face deformations. 

5. Load factor variations due to high pressures, distortions, entrance 
effects, etc. can be accommodated by incorporating self-acting lift pads to 
the seal face. This is especially true for multi- point operation. The hy- 
drodynamic force generated by the lift pads is shown to keep the seal 
"balanced” even though the secondary seal diameter is fixed. 
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Table 1 - Other names used for the geometric balance ratio or 
modulus and pressure profile load factor that are found in the 
literature 


I Geometric balance ratio or modulus - h y drostatic Pressure closing area 

sealing dam face area 

Name 

Reference 

Area ratio 

(2) 

Balance factor (used by some seal manufacturers) 

(3) 

Pressure balance ratio (modulus) 
Geometric unbalance ratio 

(4) 


n. Pressure profile load factor F = sealed Pressure force or 

sealed pressure opening force 

Pressure profile factor = net or ave seal face Pressure 

sealed pressure differential 

Name 

Reference 

Balancing modulus 

(4) 

Flow coefficient 

(5) 

Load Ratio 

(3) 

Pressure form factor 

(6) 

Dimensionless force 

(7) 


Pape (3) also defined a complement to unity which is called the degree 
of balance. 
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Table 2 - Seal force balance indicators for 
several design conditions 


Design point 
(represents aircraft 
engine operating 
mode) 

Geometric 

balance 

ratio 

Pressure profile 
load factor 

Net closing due 
due to hydrostatic 
pressure, a 
lbf 

1 - Idle 

0.70 

0.61 

4.7 

2 - Cruise 

; 


.68 

3.0 

3 - Takeoff 



.69 

2.0 

4 - Climb 

V 

/ 

00 

CD 

4.5 


cl 

Pneumatic closing force minus opening force. 
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Figure 1. - Pressure-balanced face seal. 
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Figure 2. - Schematic of radial face seal identifying the nomenclature. 
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Figure 3. - Model of the primary sealing faces (Dam) with 
a small tilt angle (not to scale). Diverging face shown. 
(Upper ring removed for clarity in top view.) 





Figure 4. - Load factor variation with pressure 
ratio in the subsonic flow regime, parallel 
sealing surfaces. 



LAMINAR TURBULENT 



PRESSURE RATIO, P 0 /P 3 cs-66446 

Figure 5. - Load factor variation with pressure ratio over a wide range of 
conditions, parallel surfaces with a mean film thickness of 0. 2 mil. 
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Figure 6. - Load factor variation with pressure ratio for subsonic flow. 
Both linearly converging and diverging sealing surfaces are shown. 
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Figure 7 . - Load factor variation with film thickness for 
several pressure ratios, parallel sealing surfaces, 
isentropic entrance flow conditions. 



LOAD 

FACTOR, 

F 



Figure 8. - Load factor variation with film thickness for several 
pressure ratios, parallel sealing surfaces, entrance flow 
represented by a 0. 6 entrance loss condition. 
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Figure 9. - Load factor variation with film thickness for several 
pressure ratios, linearly tilted sealing surfaces of positive 
one milli radian (divergent). 
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Figure 10. - Load factor variation with film thickness for several pressure ratios, 
linearly tilted sealing surfaces of negative one milli radian (convergent). 
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Figure 11. - Load factor variation with 
film thickness for positively and neg- 
atively tilted surfaces of two milli- 
radians and a pressure ratio of 10. 




Figure 12. - Equilibrium gas film thickness as deter- 
mined by total seal opening and closing forces, 
parallel surfaces (11). 



